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msruixas 


An analytioal f emulation for no derat ely large anpli- - 
tude travelling wave nenbrane flutter in inconpressible flow 
has been attenpted here using finite elenent nethod, IThe flow 
has been assumed to take place on both upper and lower surfaces! 
of the nonbrane having flexible boundaries* Ihe analysis has ■ 
been nade with the assemblage of 16 rectangular elements 
using the 12-paraneter non-conf orning model# Hermit e polynomials 

of first order have been used as interpolation functions* An | 

} 

expression for the lateral deflection, - representing the tra- [ 

i 

veiling wave behaviour for the entire membrane in the global i 
set, has been used. She inplane tension effect ap.6. the non- 
linearity due to large deflections has been considered throu^ 
the derivation of the geometric stiffness matrix. Small pertur- 
bation theory applied to the non-circulatory incompressible 
flow past the oscillating mmbrane has been used to give the 
aerodynamic loading. An aerodynamic matrix, which is kinma- 
tically consistent with the stiffness and mass matrices has 
been derived using the virtual work principle* Ihe flutter, 
criterion has been discussed and a method to solve the coE^lex 
eigenvalue problem has been presented# 
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CHAPTER 1 
IHTRODUCTIOH 


1.1 General 

During the last one decade or so, there has been 
much interest in the utilization of aerodynamic surfaces 
formed of flexible inembran(j!es. Y^hile sails are the histori- 
cal example of such surfaces, their applications are also 
found in other examples such as parawings, parachutes, fle- 
xible rotor blades, glider wings using flexible lifting 
surfaces and possibly the inflatable aerodynamic surfaces 
which have yet to come into being. The aerodynamics of the 
flexible surfaces is different from that for rigid surfaces 
since in the former deformations influence load distributions 

The aerodynamic analysis for the flexible lifting 
surfaces has received the attention of few investigators 
(1,2,3). They have established the shapes of these surfaces 
as a fimction of tension and free stream dynamic pressure to 
give the prescribed aerodynamic loadings. It has been found 
that for a particular value of inplane tension and for a 
given configuration there exists a free stream dynamic 
pressure beyond 'which the monbrane becomes unstable and' a;^ 
oscillation. becomes divergent. At this stage the flexible 
Burfh6^‘ ' ^'0'’ beOTe, ■funofioh nf- .'pr©vi4^f 
navigating force tp' the ;'.v?4hiple* Therefore, ■It ' ' ■ " 
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essential to investigate this aspect of dynamic instability, 
widely known in aeronautical literature as 'flutter*. 

In general, flutter is defined as a self excited 
or unstable oscillation arising out of the simultaneous 
action of elastic, inertia and aerodynamic forces upon a 
mass or a system of masses. During this oscillation it has 
been shown that the body extracts energy from the air stream. 
Apart from the flutter of lifting surfaces, this phenomenon 
is observed in certain non- lifting surfaces also as in case 
of flags, sails, suspension bridges, electrical transmission 
lines, skin panels forming external surface of aircraft 
and missiles etc. The flutter of skin panels, generally 
termed as 'Panel Flutter ’ is a local dynamic instability of 
a panel or group of panels of the exposed surface of the' 
vehicle and is thus distinct from the flutter of a lifting 
surface or vehicle as a whole. 


The flutter encountered in flags, sails, parav/ings 
etc. is also included in the study of panel flutter where 
the membrane acts as a panel. In these cases because of the . 
nature of boundaries and high flexibility of the struotur@r,. 
it is quite likely that flutter can occur even in imon^es** 


sible flow and with quite large amplitudes, 


observed: in. actual praotioe ■and’W-ind tunn.ei: 


-the faot‘-that, 




moving ih the th®- airs’kre^^, Ib'uS' an' aeroel^ 

ti(^^ id s . c cmfr«^ed';'isd|^i frofel^ ;''0f hi^i' a33i>lit ude ' 
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travelling wave nenbrane flutter in incompressible flow. 

There seems to bo very little literature available in 
this area, Motiva.ted by this fact, a preliminary attempt 
has been made to analyse this problem in the present work, 

btuvtM avtes 

To account for the arbitrary b aund ^ gy o ondi t-i-^'ns and to 
treat the problem as realistically , as possible, finite ele- 
ment method has been employed for the analysis. 

1 • 2 Literature Survey 

As stated earlier, very meagre literature is avail- 
able on travelling wave flutter of membranes with flexible 
boimdarios, in incompressible flow. However, the problem 
of panel flutter both in subsonic and supersonic flows has 
been investigated extensively during the past twenty years, 
both theoretically and experimentally. In the following, few 
investigations relevant to the present work will be reviewed. 

The first recognition of panel flutter as an engi- 
neering problem was by Jordan (4) who indicated that a 
number of early German ?-2 missile failures were suj^osed 
to be due to this instability. His investigations were based 
on unpublished war time research. He identified the fact 
' that a flag fluttej^ at quite low speeds forming travelling 
Waves '• and ' ao does ' 'a . skin : panel made ■ ef ■ thin' ;:Cmembna#.).;* 

if it i^ -not poQ . ti^tjy- .ifbreeh^, ClontmrF tp, the. theore- 
. tieal predictions he d fia^ of fs^ric held in, 
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a supersonic stream will no longer flutter but will stay 
stiff as a board# Shas he concluded that paper or fabric 
skins can flutter at quite low .speeds whereas .metal skins 
due to their higher modulus, require higher speeds, and if 
they flutter at all> do so with very ssaall amplitude. At 
sufficiently hi^^ speeds, any disturbance that mi^t arise 
is dissolved into smaller and smaller waves and the. result 
is a skin that gives the appearance of being stable, 

Noticing the analogy between the panel flutter and 
the classical problem of the generation of water waves by 
wind and flutter of flags and sails, Milos (5) carried out 
panel flutter analysis by studying the progressive waves in 
a plate immersed in a flow. He represented the transverse 
deflection by a travelling wave motion and solved the complex 
eigenvalue problem, Goland and Luke (6} treated the problem 
of a vibrating membrane in a supersonic stream by the method 
of Laplace transform, Hor pressure distribution they used 
the expression given by Miles (?) for an oscillating airfoil, 
and found that for the membrane, flutter cannot occur at 
hi^ supersonic speeds. Hedgepeth (8) explained the qualita- 
tively dissimilar behavior of plates and membranes at hi^ 
Mach numbers. According to this analysis, plate panels flft-, 
tter at high supersonic Mach numbers virile membranes do not# 

■ ' / 'In' ref erenc e (9) ' :ie30Pl5,ated' -;a\ 
panel' flutter problem' in the-''^'©-^ 'of '^.^.integt®'l;;'Otbation 
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and solved it nunierically by the ncthod of iteration and 
the method of matrix approximation. He also pointed out the 
difficulty in using the Galerkin method for a membrane since 
the existence of a real eigenvalue depends upon the conver- 
gence of the solution. Again Fung (10) suggested that for a 
realistic theory the panel flutter should be studied as a 
problem in nonlinea,r mechanics. 


Dugundji et al. (11) investigated the subsonic aero- 
elastic stability of a two-ilimcnsional panel reoting on a 
continuous elastic foundation both theoretically and experi- 
mentally. Ihey also observed a definite flutter of a travel- 
ling type. Instabilities were found to occur when the complex 
values of wave velocity (eigenvalue) possess negative imagi- 
nary parts. Por a finite panel Galerkin type solution was 
obtained by assuming the modes Satisfying the pinned end boun- 
daiy condition. Hedgepeth et al, (12) analysed the dynamic 
stability of an infinitely long panel on equally spaced sup- 
ports, having an airstrean of arbitrary speed on one side and 
dead air on the other; using two-dimensional linearized com- 
pressible flow theory together with elementary beam theory, 

A travelling wave expression was assumed for the transverse 
deflections. Ihe problaa of finding the aerodynamic forces 



Was reduced to the determination of the ■pressures acting qn 
a .stationary wavy; wali,.; bf,;jia^e-given,by. tbe, spatis^' jp^iq-' 
dicity of, tranSTSr'sS 
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equal to (free streani velocity + wave velocity) above and 
a flow of velocity equal to wave velocity below the wall. 

faneda (13) investigated a classical problea of 
waving motion of a flag experimentally in a low turbulence 
wind tunnel. He observed that the flag flutters in various 
oscillation modes with a travelling wave motion at the trail- 
ing edge. Sunderarajan (14) studied the effect of a rigid 
boundary placed parallel to a long panel with finite width 
assuming a travelling wave expression for transverse defl- 
ection. The aerodynamic pressure was calculated using Courier 
transform technique. He concluded that the effect of the 
rigid boundary is considerable in lowering the flutter speed 
when it is placed very close to the fluttering panel. In ref- 
erences (15)> (16) and (17) the supersonic membrane flutter 
problem is considered as a special case of a panel with neg- 
ligible flexural rigidity but having hi^ inplane tension. 
Dowell (18) and Morino (I9) studied the problem of nonlinear 
panel flutter. The nonlinearity due to the membrane forces 
was included in Ton Karma’s large deflection equations. 
Dowell applied Gelerkin method and the system of nonlinear 
ordinary differential equations in time was solved fey numer- 
ical integration^^. 

Hie method of analysis in the above Investig^ldLons 
fetK>ad;^, .coijsists of .the f oll,OWing 
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(i) Expand the structural deformation in a series of 
natural or atleast complete modes. 

(ii) Determine the aerodynamic forces for the given 
mo d al d ef o rmat i’.-) ns . 

(iii) Using these in the equations of motion, apply Gale- 
rkin’s method to arrive at a set of ordinary, nonlinear, 
integro-diff erential equation in time for the modal amplitu- 
des. 

(iv) Solve these equations by a numerical integration 
with respect to time. 

In the literature there are many individual variations 
from investigator to investigator in the use of theories to 
get aerodynamic loading, in the representation of the panel 
deflection and deriving the equation of motion. Is structures 
became more complicated with arbitrary boundaries, the con- 
ventional methods like Galerkin^s and Rayleigh-Ritz proved 
incapable of yielding accurate solutions to them. Recently, 
therefore, few investigatCJrs ( 20 , 21 , 22 , 25 , 24 ) have attempted 
to use finite element techniques to solve flutter pioblems. 

Eariappa et al. (21) analysed successfully a simply 
supported skew panel in supersonic flow tlft^gh, this approach 
Aerodynamic influence coefficient matrices that are kine- 

#ere, derived usi3Q.g Of virtual werki use 
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of assumed modes was avoided by solving? directly the complex 
eigenvalue problms for studying; tiie stability characteri- 
stics. A further application of finite elements was prese- 
nted by Olson (22) in which he used two rectpjigular plate 
bending elements (l2-paremeter non-conforming model and the 
16-parameter conforming one) and an 18-parameter conforming 
triangular element. He concluded that conforming rectangular 
element yielded higher , aocurany and faster convergence than 
the non- conforming one. A further extension of work of (21) 
is presented in (23,24). In (24) Kariappa et al. considered 
the flutter problem of skew panels with inplane forces under 
yawed supersonic flow, (Dhe effect of mid plane forces was 
included using the concept of geometric stiffness matrix. 


1 *3 Object 


Present work deals with the analysis of^membrane flutter 
of travelling wave type in an incompressible flow using fini-. 


te element technique. The analysis made is completely a g<aie- 
ral one and is valid for any set of kinematic boundary condi- 
tions. The main ob;}ect of the present work is to provide a 
theoretical analysis for the aforementioned problem, which 
has been studied experimentally in the Dept, of Aero, Bngg, 
I.I.T, Kanpiic* The idea being to predict the results 
tically. A pheok on the analysis Pan be made pnd n®Gessar;y \ 
mo^fications be inco3:poratsd f fb needed,'. by :'pQmparing bhpse ' 
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theoretical results with those of experinental study. In 
the tests the rectan£;ulRr convasSi cloth model was tested 
by fixinj:? it at its four corners. Those boundary conditions 
were chosen, because the fluttei* could occur only for these 
conditions within the available tunnel speed range. Therefore, 
in the present analysis also this set of boundary conditions 's 
considered. 

Chapter 2 gives an account of ejxpcrinental . study of 
the problem. It includes the experimental set-up, instruments 
and procedure. Observations and inferences have been presen- 
ted t.i act as a guideline for the theoretical analysis. 

Chapter 3 gives the structural idealization, of the membrane 
by dividing it into 16 rectangular elements. The function is 
discretized using Hermite polynomials as the interpolation 
functions. Kinematically consistent aerodynamic matrix is 
derived through virtual work principle by representing the 
menbrane with a sheet of point sources and using linearized 
small perturbation theory. In Chapter 4 the stiffness and mass 
matrices are derived using the principle of strain energy 
and virtual work concept. lonlinearity due to hi^ amplitudes 
of the fluttering membrane and the effect of inplane tension, 

* 

are takeaa into account throu^ the derivation of geometric 
stiffness matrix. Transformation matrices from local to gplobal ' 
cooMinsft; es are derived: :in. :C&apt ^ ■■5 ^and method. is;''-^yeil,^;';\-;c 

for '^sembling the ';aQro^dyhami.-C^,; stiffness, .and 'mass, matri:e'^;.#t::;:>: 
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iifter accounting for the boundary conditions the complex 
eigenvalue problem is fornulated* Chapter 6 indicates the 
method of solution. Appendices at the end gives the detailed 
calculations of stiffness and mass matrices and the trans- 


formation matrices* 



GH.^TER 2 


EXPERIMEN-Iia Sl'UDY OF THE MEIffiRaiME 
EHJIgER PHENOMOII 

2.1 General 

As described in Section 1.1, a necessity was felt to 
study the phenomenon of nembrano flutter. A theoretical ana- 
lysis Was therefore needed to predict this behaviour. To 
make this analysis more realistic, it was decided to have a 
physical picture of the flutter phenomenon, so that it could 
act as a guideline for the theoretical investigation. Simple 
experiments v? ere therefore conducted using the tunnel avail- 
able in the Dept, of Aero. Engg., I. I.T. Kanpur, The necessary 
apparatus was built by the author and Santa, who was interes-' 
ted in making very detailed parametric studies of the pheno- 
menon (26). 


2.2 



2.2,1 Wind Tunnel 


The experiments were conducted in the low, speed 3*x2*, 
three dimensional closed circuit wind tunnel, available in 
the Department of Aeronautical Engineering, I.I.T, Kar5)ur, 


The maximum speed of the tunnel was 42 metres per sec- 
ond. To observe the. amplitude and shape of the flat terlrsg 
membrane a grid work; of 2 ems sfuare was, made on a slieet of 
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black paper and was pasted on the facing- vertical v.'all of 
the 5 ft, long: test section. This also provided the necessary 
black background for the model to take the photographs. 

2.2.2 Design of Membrane Fixture 

While designing the fixture, due considerations were 
given so that nininun alterations were required in the exi- 
sting test section. The other factors considered were the 
ready installation and ease in fabrication of the fixture. 
With these considera.tions a very simple and multipurpose 
fixture was desi^ipaed and is shown in figure (2.1). If consis- 
ted of four adjustable, each 60 ems, long and 2 cns. in dia- 
meter, vertical mild steel rods threaded all along their 
length, with 25 cns. long flanges attached at their top. In 
the test section floor six holes, of 2 cns. dia. each, were 
drilled. The vertical supports could be fixed in any set of 
four holes corresponding to the required dimensions of the 
membrane. Additional extension strips of various sizes, as 
shown in Pigure 2.i, were provided to accommodate any confi- 
guration any dimension and any type of boundary condition* 

2.2.3 Test Model 

S'or observation purposes a rectangular membrane model 
was prepared out of a 0.06 cms. thick convasi| elolii ^eet^ 

' '--In i actual pract±'ce.-f he _,reQui3?^ents;of -edge , hondi.tions , 
for the membrane ahe ;diffe'a?eht .'^plication^ 
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was possible to test the moahrane with almost all possible 
type of edge conditions with fixture designed in Section 
(2.2*2). Trials wore performed at the following sets of 
boundary conditions. 

(i) iill edges clamped. 

(ii) Only leading edge and trailing edge clamped 
and sides free. 

(iii) leading edge clamped and trailing edge fixed 
with flexible mild steel wire, sides free. 

(iv) Both loading edge and trailing edge fixed with 
flexible mild steel wire. 


In all these cases the membrane vjas stable within the 
available speed range. 


Then it was decided to further relax the boundary 
conditions by fixing it at its four corners only. In this 
case all the edges were stiffened by applying a thin coat 
of dope on them, ^‘or this set of flexible boundaries the 
flutter occur ed within the tunnel speed range. So, these 
conditions were retained for further observations. 


2.2.4 Pressure Probe 

It Was planned to record the pressure distribution 
over the membrahe durii^ flutter with the idea of determi- 


ning the pressure distribution in the nei^bourhood and at 
instability.' Por this 'purpose a pressure probe, was made 



which consisted of seven static and one total head tubes, 
covering the entire v/idth of the nenbrane. iPhese tubes were 
attached at the ends of two long brass tubes st arming from 
the top surface of the test section (as shown in Figure 2#SI), 
:%■ traversing that probe from leading edge to trailing edge, 
the pressure distribution on the complete membrane, both 
along stream and across stream could be recorded* 

For this purpose a traversing arrangement, as shown 
in Figure 2.S, was devised, lo fit this arrangement the exis~ 
ting top of the tunnel test section was replaced by a new one 
having two slits of one cm. width each and 45 ems, apart, Do 
enable the traversing on either side of the membrane surface, 
while the tunnel v/as running, the probe was made in two parts 
One contained 3 static tubes and the other had 4 static and 
one total head tube. Both parts could be engaged and dis~ 
engaged as desired with the help of a *0’ clamp provided for 
the purpose, without stopping the tunnel, 

2.2,5 Manometer 

An inclined multitube manometer with distilled water 
as its working medium, was used to measure speed and pressure 
distributbn, Eunnel static and total pressure connections 
were given to the manometer so, that their difference gave the 
dynamic head inside the tunnel. Measuring the temperature 
of 'the. airstream inside, with an alcohol thermometer, the 
air densi^ corxections were made while oal^iating t^ 
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airspeed from the noted dynamic pressiire. 

The pressure probe tubes were also connected to the 
same manometer and from the nanometer readings the pressure 
coefficient was plotted along and across the membrane length, 

2.2,6 Stroboscope 

The frequency of the fluttering membrane was obtained 
by using the strobotac, y/ith this the oscillating manbrano 
was made stationary and the node shape, amplitude and the 
nature of the oscillations were observed by taking a photograph 
of the stationary menbrane as ^own in Figure 2,4. 

2.3 Observations 

2,3.1 Development of the Flutter. Phenomenon 

The membrane was fixed at its four corners on the fix- 
ture and the tunnel speed was raised gradually from zero 
onwards. The nature of the development of oscillations was 
noted at successive increments in speed. These observations 
are described in the following. 

At very low speeds the menbrane stayed in a sli^tly 
deflected shape with occasional gentle oscillations of small 
amplitude. As the speed, was raised, these oscillations grew 
larger in aap^tude and, at times j irregular sudden jerks 
were observed* At , still hi^en spf e^, m^brane took a shape 
forming .■';a;pe.ak n0W..;tlae -leading^r a fallow and long 

;dppese,iQn i^J'ailing ' ■ same speed they 
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interchanged their positions by occasional up and down move- 
ment. A slight further increase in the airspeed pushed the 
leading edge peak to the downstream side and this resulted 
in a quite regular motion of the membrane. In these regular 
oscillations j the peak formed at the leading edge, due to 
some flow disturbance, appeared to move towards the trailing 
edge. The whole phenomenon looked as if some small spherical 
caps were being conveyed from leading edge towards trailing 
edge tlxrou^ the membrane. These caps were the disturbances 
created at the upstream side, that were moving towards the 
downstream edge of the membrane. This was identified as the 
’travelling wave flutter’. The stage at which this phenomenon 
oc cured Was defined as the ‘onset’ of flutter. 

Once the flutter occured, the strobotac was s^'y itched 
on and the strobotron lamp .was focused on the oscillating 
membrane. In the beginning, the strobotron li^t was projec- 
ted at its hipest frequency. The frequenqy was then dec- 
reased gradually till the first single stationary image of 
the fluttering membrane was obtained. To differentiate betwee 
single and multiple images a coloixred spot was marked on 
the membrane surface. The single image was observed when only 
one image of the spot was seeio stationary. While making the 
meabrane stationary the waves appeared .moving against the 
at|r.-e^ , -at the strobotac. ^settiiJg. higher' '-than ,thq ' ; 

'queucy. At a lower settin^,;^,#y.'^p,eared;ja5.Tlti:->al.qng'-the-:,, - 

''' "'"V , " 



17 


stream. In between these two readings there was a setting 
at which the nsiibrane jnst appeared to be stationary as 
diovm in Figure 2*4* At this condition the node shape consi- 
sted of generally 3 peaks and depression (in some cases there 
were 2 peaks also). Ihis setting corresponded to the flutter 
frequency and the corresponding tunnel speed was the flutter 
speed. 

2*3 Pressure Measurement 

Knowing the flutter onset from the first run, the 
airspeed in the second run was kept to a value slightly less 
than the onset speed, fhe pressure distribution was measured 
at this speed. In the manometer, the pressure probe readings 
showed a decreas.e in pressure over the peak and a relative 
increase in pressure over the depression on both the surfaces. 
Next, the airspeed was raised till the flutter speed was atta- 
ined. Again the pressure distribution was measured at this 
stage. Ihe pressure plots are shown in Figure 2.6. 

Because of the unsteady oscillations, the readings 
in the- manometer were also changing from instant to instant. 

At any instant the. pressure plot was resembling veiy nearly 
the deflected shape of the membrane. As the membrane m^de 
shape Was changing, so did the pressure distribution, lo 
record a time history pf these unsteady pres si^e reading : ; 

hi^ speed-.; f'ilm, records, W;^S taken with- a mp^vie.'casenai 

;But the d.ey.®;lrOp®i^t pf;;thf:rf.i3li^oanid,-^b3t'' be .arranged .in, time:. 
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So, only the mean manometer readings were taJcen and corres- 
ponding mean pressure distribution was plotted to get a 
rou^ idea of the a,erodynanic pressure distribution over 
the fluttering membrane. 

2.3.5 Inplene Tension Effects 

Some observations were made to see the effect of 
Inplane tension on flutter. Por this the membrane was stre- 
tched arbitrarily in the beginning and then tested. In that 
case it was observed that the flutter speed was raised and 
the lateral deflections were also restricted, A relatively 
lov/ amplitude travelling wave flutter was observed with 
hi^er flutter cpood and frequcncs’’, as shown in Pigu-re 2.5. 


2.4 Inferences 

From the observations made in Section (2,3) following 
inferences were drawn* 

1. luring instability, the membrane exhibited waves that 
travelled towards downstream side. Ihus the membrane flutter 
phenomenon Was of travelling wave type. 

2, Edge conditions played an important role in the membrane 
flutter behaviour, Ihe more the edges are restr^ned,. the 
higher is tide flutter speed,, 

3* ■" la’;. ee^aih^G!a^ses: : membrane .Was flattening 

others, quite l^ge , amplittudes;' , 

4V-' Inplana ;t^pibia;-ha#;,;,a’:fayo^ effect-; ''in the , 



CHiiPTER 3 


STRUCgURiJ. IDEALIZilTION OP MEMBRi^KS im 
DERIYaTION of ilERODYRMIC MiiTRIJL 


3.1 Einite Element Method 

In any elastic continuum the solutions for both the 
static and dynanic prublems, arc in most cases very difficult 
if not impossible, because of infinite degrees of freedom* 

!Ehe concept of finite elements attempts to overcome this difi- 
culty by approxima,ting the real continuum to a finite degrees 
of freedom. In this case it is visualised as an assemblage of 
simple stiuctural elements interconnected at a, discrete num- 
ber of nodal points at which some fictitious forces are sup- 
posed to act. Ihe problem then reduces to the so'lution of a 
set of these simple element s» amenable to numerical analysis. 


This concept of finite elements can briefly be sum- 
marized in the follov;ing manner* 

(1) The continuum is divided by imaginary lines or 
surfaces into a number of finite elements. (Discretization 
of Region). 


(2) The elements are assumed to be interconnected 
at a discrete number of nodal points situated on idaeir 
daries. The displacements (or the forces) at these h|o,del poihti ;! 


■will be the .b&sic,unkiKjwn ;par#i#ers’ of .Itie.p 


robies* 



20 


(3) A function is chosen to define uniquely the 
state of displacanent (or forces) within each finite element 
in terns of its nodal displacements. (Discretization of 
function) 

(4) Ihe displacement functions now define uniquely 
the state of strain within an element in terms of the nodal 
displacements# These strains,, together with any initial 
strains and the elastic properties of the material .will def- 
ine the state of stress throughout the element and, hence, 
also on its boundaries. 

(5) A system of forces concentrated at the nodes and 
equilibrating the boundary stresses and any distributed loads 
is determined, resulting in a stiffness matrix for the element 

(6) Knowing the stiffness properties of individual 
elements, the behaviour of the assembled structure can be., 
studied throu^ any well known technique of structural ana- 
lysis. 

Thus finite element analysis may be thou^t of as an 
assumed mode analysis in which the assumed modes, instead 
of having an uniform definition over the entire structure, 
have a different definition in each zone or subelement, 
producing a sort of discretised assumed mode. The local modes 
(usually polynomials) ,;are assumed, for convenience, su<h 
that"ti|^^.Ulldeter^ned,c;o'^f|,eieh|e';.edhhe^.ond. to; physical 
d-egrees of ^eedom o,£ elaaent,# ^ ■ 
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!Elic finite elonent approximation essentially reduces 
to the pnoblcn of miniriizing a- *functional* representing 
the total potential energy defined in terms of a finite 
number of nodal pa,rameters and formulating* a set ,pf simul- 
taneous equations 4 iimong all the available techniques of 
structural analysis, finite element approach is considered 
to be the most powerful technique and is best suited to ass- 
oable the complex and realistic configurations from relatively 
simple element shapes. 


3.2 Structural Idealization of Membrane 
3.2.1 Biscretization of Region 

The degree of aj^proximation achieved througji the 
finite element method depends very much on the choice of 
element shape and of the form of the displacem« 2 it function, 
The choice of element shape further depends on the nature of 
the problem. In vibration and flutter problems of panels the 
plate bending elements are used. 


To get a real picture of the problem this choice was 
based on the experimental observations described in Chapter 2, 
These observations can be of a great help in selecting the 
shape of elements to approximate the exact deformed confi- 
guration of the membrsnci 


3.2«1>1 C^Qic.e:-Qf Memenf'i^aape '.' ' 

figurations, '' f lat ■ 't ^quadrilateral ,ei;®'e3its give 
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the "best approxination. But in the case of regular confi- 
guration like a rectangular panel, the simple flat rectan- 
gular sha,pe is the obvious choice, 

Bor oscillations of moderately large amplitudes as 
shown in Figure the deformation of the membrane remains 

within the reach of shallow shell theory. Thus while applying 
the incroaental step approach, the behaviour of this conti- 
nuously curved surface, formed after first step, can be ade- 
quately represented, during any subsequent step, by the be- 
haviour of a surface built up of small flat elements. The 
only difference in that case will be that local coordinates 
for these flat elements will keep on changing from step to- 
step. 

Bor still larger deflections as shown in Figure ( 2 ^ 
the ^ape of the membrane becomes more shell like and the 
above representation seems to be unsatisfactory for this case. 
Here the curved shell elements have to be considered in each 
step to approximate the deformed shape of the membrane, 

3. 2. 1.2 Oho ice of Element Size 

The investigators in the past have shown that in thin 
plate structures even a gridwork of 4 elements (2x2) gives 
reasonably good results. As the size of the element is de-; 
creased, ;the uceuraey, of course, increases but the mathe- 

'' iiiTOj,y®d.:" As' utSserved,' in' these ■: figures : ■■■': 
'(2,^ .'l.6,elments; ' 



forming a 4x4 mesh is expected to give good results. 


3. 2. 1.3 Division of Membrane into Discrete Elements 

In accordance with the considerations discussed in 
Sections (3.2.1 .1) and (3 .2. 1.2), the rectangular membrane 
shown in Figure 3.1 is divided into 16 flat rectangular 
elements. Thus a relatively finer 4x4 gridwork is obtained 
as shown in Figure 3.2. The local co-ordinate axes for a 
typical element, as shown in Figure 3.3, are taken along the 
global co-ordinate axes. With this selection of co-ordinate 
axes the elemental characteristics obtained in the local co- 
ordinates are the same as in the global co-ordinate system. 
This helps in simplif^dag the assembly process. 


3.2.2 Discretiza,tion of Function (Choice of 

Element Shape Function) 

The state of deformation of the oscillating membrane 
can be described by the lateral field displacement w, which 
has to be expressed in terms of nodal displacements using 
the element shape functions. At any node 'i* of the element 
the following generalized displacements are considered. 

Lateral displacements and out of plane rotations 
i.e, dw/8x, Qw/^y. 


Suitable shape functions have to be detemini^ 

Two types; of sh^e'fu^ are deiihed.as ibllo'^s?,. ' 

{if. '■ ' sb^e ; fnnetidiis. wHieh "'eiMiljes.:: ihe coidpiefe'' 

continuity of di^piaoemi^t w of siibpes 9^/9 x, 9w/9y sLong 
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an interface between various eloaents. 

(ii) Non-confonaing shape functions s which preserve the 
continiiity of w but violates the slope continuity between 
the elements, though not at the node where continuity is 
imposed. 

In the solutions of first type the mathematical and 
computational difficulties often rise disproportiomtely fast. 
However, it is relatively simpler to obtain the solutions of 
second type. 

Proceeding from the idea that the imposition of slope 
continuity at nodes must, in limit, lead to a complete slope 
continuity, several very successful elements have been deve- 
loped (Zienkiewicz, Glou^, Argyris etc.)* Bazeley et al. (27) 
have established the superiority of non-conforming solutions 
over conforming ones for most of the static problems. To verify 
this superiority for the menbrane flutter problem analyses of 
both kinds have to be made, m the present analysis, the non- 
conforming shape functions have been used. 

By using the products of one-dimensional interpolation 
functions it has been observed that it is possible to goie- 
rate displacanent states having the property that displace- 
ments along an ed^ depend only on degrees of freedom c&mmn. 
to/that, _edge.''2his 'is';;all,-tha‘fe ls,_ req.ul3^a 
oompati^ilty bet^'.^enj'the' v- 
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functions selected here are the Heraite polynomials as used 

/ 

by Bogner et al.(28). Retaining idle same notation, a Hermit e 

polynomial Hj, (x) is a polynomial of order 2Hf2 which gives 
i 

when X = x. 

d“H/^dx^^ =1 m - k for k - 0 tc H 
and d^/dx^ =0 m k 

where H is the number of derivatives that the set can inter- 

polate and x. are the specific values of the argument x of 
J 

the polynomial. 'They can be used to interpolate a function 
f (x) and its derivatives at the two points *0* and given 

the values f(0), f(a), f*(0) and f^(a), as follows: 

f(x) = f(0) + f(a) + fi(0) ^ f 

In the present case f(x) corresponds to w and f(0), f(a) 
etc, are the Nodal displacements Yfj^s ett. In general they are 
defined as 

S H(p(x) = 

13 

2 

H^^Cx) = (x-x^) [h^^^(x)] .... and so on. 

Ror the present case of flat rectangular element with sides 
a and b the polynomials required are defined as follows: 

H^l^(x) <=: ^ ( 2x^-3 j « •^2y^--3by^4-b^) (3.1a) 

=: -,h (2s?-r3as^) j Edi5(y) - - XrXzp-'i'^r\ 0.ib) 

^ C3.''»1c) 





26 


H^^Cx) = ^ (x^-ax^) ; ^ (y^-by^) C3.1d) 


Vfitli the help of these polynomials the following' 
field displacement is generated. 

Assuming a travelling wave expression, the transverse def- 
lection is written as 


w 


i2VA (ct-X) 2 2 

(x,y,t) = e I Z 

4—4 4—4 Ox 


i=1 3=1 


(x) ®i3+H^?^Cx), 


V ) (3.2) 

Where A s= v/ave length 

c « wave velocity 

X is the global x-co ordinate indicating that the travel- 
ling Wave phenomenon occars in the entire assembled membrane 
and not in each discrete element separately. 

She inplane displacements accompanying these lateral 
oscillations can be neglected, because their effect being 
of second order, this does not cause any substantial error. 
Where 'W. . etc, are the values of w at the point x = x. j 
y = y. i.e. w . = w(x. , y.) etc, (Refer to Figure 3.3) 

3 x j 

In the matrix notation Eqn, (3.2) is written as. , 
follows* „ , 

w(x,y,'t) = e . ■(^! . (3.3) 

where A; is the wave .number by X - 2Vx C3«4) 
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[w^ Wg W3 W5 Wg Wg Wg w^g] 

(3.5a) 


h(i) 

^01 

.(1) 
; 02 

12 


(x) 

* 

HoV<y) 

(x) 

m 


(x) 

m 


(x) 

0 


(x) 

# 


(x) 

• 

H^^^Cy) 

(x) 

« 

H^l^Cy) 

(x) 

• 

H(J^(y) 

(x) 

• 

H:[|^(y) 

(x) 

• 

HQ2^(y) 

(x) 

• 

H^>(y) 

(x) 

• 

H|^^(y) 


(3.5b) 


5.3 derivation of Aerodynamio Matrix 

3*3.1 Press'ure' lietribution on the Membrane Element H 0 .I 
Tile pressTrce distribution on an oscillating .iftng or 
airfoil subaerged'/in,:' unst^^j OQsaxrossible ■flbW'‘>is’-'Obt''i^'4oi'' 
by, solving^tfee iiao^ife#€t;'i^a^iai:;^ffe2?''^tiai''''Oiiuatib3i-"ife 
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disturbance velocity potential ^ . (Refer to Bisplin^off 
et • al ( 29 } • ^ 



(3.6) 


Subject to the linearized boundary conditions 

§z ” Y 5 for z = 0”^ (x,y) in (3.7) 

fk~ ^ lx * z 5= 0 (x,y) in (3.8) 

where Ra is the portion of the x,y plane covered by . the 
projection of the planform. Biua'|;ions (3.7) and (3.8) repre- 
sent the exact specification of tangency of the flow to the 
surface. For incompressible flow, the speed of sound becomes 
relatively large and equation (3.6) reduces to the Laplace’s 
equation 

^ 0 ( 3 . 9 ) 

This method has its foundation in the assumption of small 
disturbances. Quasi steady state assumptions ^e made which 
amounts to say that entire field of particle velocities de- 
pends only on the instantaneous motion of the submerged mem- , 
brane, qtdte independently of the past history of that motion. 
!Ehis behaviour stCTs from the fact that the speed of sound 
is effectively infinite, so that any chgtnge of bonnidairy con- 
ditions is propagated; fnstantan'equsiy to.'' all ^^i^i.dlesf 
Here 'a 'nop5ircmlatO:xy.;;iinc,€i^eesibl,e'' pet-'shtial' -consi*- 

deredf ^in '"the^’above ■■ 
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mentioned property and is easier to deal with. 

Thus our prohlcis now is to solve the partial differ- 
ential equation (3.9) subject to the boundary conditions 


8l 

« 91+ 7 ^ 

6t ^ 

(3.10a) 

on R 


all 

iz=o- 

a 


« (0|. + 7 

^at ^ sx*' 

(3.10b) 

|1 = 0 
az 

off R 

B. 

(3.10c) 


where is the region consisting of the monbrane element 
placed inxyplene and w(x,y,t) is the transverse def- 
lection of the membrane. The porosity effects are .considered 
here. 


familiarity with the elementary tools of fluid mechanics 
suggests that all these requirements are met by distributing 
over R_, a sheet of those solutions of Laplace *s equation 
known as point sources. A single -concentrated source centred 
at X , y ,2=5 has the velocity potential 

- H 

- " -■- • ■ -irrrr - — ,;.n (3.11) 

4iCA/(x-f)^+ (y-"n)^+ <z-i 

Equation (3.11) describes motion with spherical ^nmetry 
around its centre, from which 'the source is readily ^osm \ ' 

to be emanating a bulk: of liquid into the surrounding space 
' equal ''to V 'Its ■ strength; Hi A/'i^eet of these,' spicead 'eo’ntimoui^^' 

, over' the: f .lirface' ;SC’ 'f pen ■' unit ^ea- 
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in the neighbourhood of the point ( 0) possesses the 

disturbance potential 


Now 


fi(x,y,z) ^ - ~ I f - 

a 

to satisfy tho boundary 


H(tn) 4 d»^ 

- — 1 — ■ (3.12) 

■>/(x- ? )^ + (y-"^ 

condit io n (3.10), we c alculat e 


aje I 

9 z • 2=0 


1_“f# ^ r r ' ) 'll a>i 

4Tt Hzi-' r 2 2 r 

^ vCx-f r + (y->) + z^ 


lahing tho differentiation inside the integral ^ign, 


as limits of do not depend on z, we get 


H lz»0+ = 4it 


1 aw. . , H(t.-n ) ar a>i 
^ 2-<rz / / 


R, 


a 


x- 1 ) + (y-^ ) + z ] 


3/2 


(3.13) 


As z assumes smaller positive values, the integral is caused 
to vanish by its multiplying factor, except in the vicinity 
of point ^ x,^= y, where the integrand tends to infinity. 

iEiis region is isolated with a small square of side 26, obtain- 
ing 


\ z * z*=0 


7 ^- z-Cf- S S 

y-e x-e 


H(f,->1 ) at dll : 


(3.14) 


L(-f)"+(y-v,)^+.f^ 

Being a continuous function, over the entire square 

differs from its central value H(x,y) by an amount of order 6# ; , 
Hence, if we neglect these small variations and intro.duce 
tQBporary'.-lnfc.egratiqn./'vari^les'.C*- ,Cr-^5f 
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Si 1 + 

dZ '21=0 


H(x,y) 

— _ z . 


? / t* ^ 

’'o rJ2 *2 ”2 tV2 
-e -e j_ +z ] 


/ 


Using the integral formula,e, 

-I- X 


dx 


and / 


dx 




o* 


'+c'xya+cx^ 


1 

(MMHI 

ci 


■/ 


T 

r 


a 


I t„ 

ac -a c 


—1 

tan X 


ac --a c 

. (a+ox5 


the above eq.uation becomes 


al., 

A <7. I 


9 2 ; * 


H(x,y) ^I»t.^ ^+ 


A-% 


2tan’'^ ( ^ 


2 ^2 
. — "I , —6 

J-2tan ( 


z/ae^+Z 


) 


;^2e2+22j 
( 3 . 15 ) 


Now z and 6 should approach zero in such a way that 
the ratio S/z becomes indefinitely large. To ensure that 
terms omitted in going from (5.14) to (3*15) are truly negli- 
gible, we have to mfike 6 vanish in such a way that v;g don’t 
get any non-zero contributions from the Region outside the 
small square. This is possible when 6 remains large comp,ared 
to z. In that case the argument of inverse tangent approaches, 
infinity because of z being a small quantity as compared to 6, 
the denominator has a strong tendency to approach zero, as 
both e and z vanish. Thus limiting in this way, the inverse 
tangents approach ^ ^ and - j' • 
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Coming up to the source sheet from the lower side reverses 
the algebraic signs of the arguments of the inverse tangents 
in (5.15) thus leading to 


9z ‘ z“0 


H(x,y) 


(3.17) 


Therefore, it is seen that the source produce the symmetrical 
discontinuity specified by equations (3.10a) and (3.10b) 
which can both be satisfied by setting 

H(x,y) = If - H 1^^- 

. ^ V g l Mi . l l) (3.18) 


Substituting (5.18), in (3.12), we get 

. a .b ajiii 


1 ^ f* 

(|)(x,y,z) = - ^ / 


0 0 


//(x-f )^+(y->} 


— df dv| (3.19) 


Pressure distribution is. now calculated by using the unsteady 


linearized Bernoulli *s equations 

AP = Pp - Pj. = -2f [||+ 7 ||] 


(3.20) 


Differentiating (3.19) with respect to x and t and setting 
z*so we get 

. .9w(t,n,t) aw(T.>i, t) 

•AP * f (|^ / / ' ' ^ 

° ° /(^-f )\(y->) )^ 

. ^ : 0.21) 

^ i i 
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Now for element no. 1 Figure (3*2) X 
(3.3) takes the following form 

i A*(ct-x) f .T 


w(x,y,t) « e 


[t>\ [f| 


X and • , equation 


(3.22) 


The aerodynamic lift force acting on the small area dx dy 
(Figure 3.5) is 


= Ap <3x dy 


(3.23) 


Substituting f or w from (3.22) in (3.21) and then calculating 
AP from (3.23), we obtain 


.1 (l2Tz/s)ot f, ,,±C>Ttc^ 2 « , 

&P - e -I- (( 1 ^) y)] 


+ (i^T) , 


+^j(2,y3|, ) + 7^ 

(3.24) 


Where the vector ‘s are given as followss 

= / / — an 




(3.25a} 






^ ^ ^ (8f/dt) - i- A* [f] 

^ M ^ ’ 


Ar^_£fl^!LkL 


dt d>| (3.25c) 






:C: ib;. 








(5.25d) 
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3 • 3 . 2 Derivation of lerodynaaic Matrix ; 

To consider tlie effect of this aerodynamic force AD, 
the aerod3momic matrix is calculated as follows: 

The virtual work of the aerodynamic force AP is 

V = /"” / w AP^ dx dy (3.26) 

^00 

where 'w is the virtual displacement. 

Now ' 

^ f^>T . , i A (ct-x) 

^ = \lf\ ^f| e 

Substituting this w and Eqn. (3.24) in Eqn. (3.26) with 
X as = 2ii/a » 'ws get 


a 


8* 




gi2 A*ct 

^ o' 

/ Wf 

0 


X*o)^ 


+ (i A* C 7) ( 

pal ' - 

I 

) 


dx dy 





(3.27) 

can be written i 

as 


'* 


> = ® M W 



(3.28) 

i the aerodynamic 

' 

3 matrix 

2 

[aJ is given by 


12 Act f a 
®* ® r / 

b * 2 

S (~ X 0 

0 

-i XX 
e 



+ i Y e"^ 

<w* 

z ( 

T 

1 P2! ^ 



^ ^2 A^x 

• f Y e ■ 


)''\'dx'"dy;‘' 


(3*2^) 


; .f'" pf,..: ^ ' dk-' *'■ i 
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This can further he v/ritten as 


fA^j = - N 0 j-A^^ 

where 

N 
N 


I + i Kq Y fk 


+ Y‘ 




(3 •30) 


i2 Kct p a b -i A*x. 
e f f e 


T 


% j j ^ \t] llJ clx dy (3.31a) 

0 0 ' 

12 Xct f ah -i A x 
e “ir / / e 


o o 








e 


12 Act P a "b *-»i A x 
r r r 


TC 


/ / e 


o 0 


^ {^2\ + (^3! t^-y 

(3.31b) 

■{f] ^ dx dy (3.31c) 


Equation (3.30) gives the aerodsmanic natrix for element no,1. 
3.3.3 Aerodynamic Matrices for other Elements 


The global streanwise co-ordinate is related to the 
other elements as follows: 

Por element nos. 1 , 5, 9, 13 j X « x 

Por element noe. 2 , 6, 10, 14 ; X »= x+a 

Por elegient nos. 3, 7j.11,15 ; Xs= x+2a 

Por element nos. 4» 8, 12, 16 ; X *= z+3a 


(3.32) 


IThe expression for w will be xcodified .according to P^n,(3.32) 
and so will the aerodynamic matrix be. The aerodynamic 
matrices for rest of the elements are p elated to that for . 
the 1st element in the foil o^ng manner: 
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W' “ n' = N 

n' - n' = N 

p^o]' = n' = M 

r -1^ 

where lA^j denotes the 


13 


1 

SC 

N 


10 


.14 -i2 Xa 

tx j 

I 


11 


1 5 ~i4 X a 

SS 1 

.■^ej 

* c 

12 


1 6 — i 6 X 3. 


K 

0 

= e f 

aerodynamic matrix for 


N 


W 


(3.33) 


element . 



CSHiJ>I!ER 4 

IlERIY Aglow Qg SglEPEBSS jA'D MilSS ILiTRIOES 

4.1 General 

In the (lynajtnio problem of elasticity, a structural 
system is in equilibrium under the action of internal forces 
namely elastic, damping and inertia forces and the external 
applied forces. In the present problem, we assume no damping. 
The aerodsmamic influence co-efficient matrix representing 
the external loading has been derived in Chapter 3. To repr- 
esent the elastic and inertia forces the , stiffness and mass 
matrices v^ill be derived in this chapter. 

4.2 Derivation of Stiffness Matrix 

In deriving the stiffness matrix, geometric non-lineax'- 
ity caused by the largo deflection of the oscillating membrane, 
is considered. This is connected with the strain-displacement 
equations. Even if strains remain small in the conventional 
sense, rotation of the element adds nonlinear terms to the 
strain-displacement equation. A geometric stiffness matrix 
is derived throu#. the use of incremental step procedure to 
take this nonlinearity into account. 

In general the stiffness matrix .can be derived using 
any of the following methodss 

(i): \nnit .■diapiaobm«^t\,: or virtti^ work' method;, ■ 
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(iii) Differential equation method. 

(iv) Inversion of displacement force relationship. 

In the present work Castigliano • s first theorem is 

used in v/hich the strain energy U is written in terms of 

f ? 

nodal displacements |W| , Ihen by this theorem the stiffness 
CO -efficient is obtained as 


k 




13 3W^ 


(4.1) 


V(/here is the column vector of elemental nodal . displa- 

cements gr.lvon by Eqni (3»5a). 2he equation (4*1) is written 
in a more concise form e,s follows s 

I PI M)f} (4.2) 

!I]ho stiffness matrix fkj , whose elements are 
is then given by the triple matrix product as 


T 


M = [A] PI W 

Now strain energy is written as 

n.. i ; ^e| ® ^<r|av 


(4.3) 


(4.4) 


where the stress and strain vectors are^ defined as 

Assiuning the material of membrane to be elastic, the consti- 
tutive law. is writtenr-S-s 


cr 

X 

J 






LI 






59 


where v is the Poisson*s ratio. 


, . For element no. 1, Eqn. (4.4) now hecomes 


i "PVi 3- b , 1 1 5) f ) 

■q / S jel le| dx dy (4. 6a) 

Vll) 


where h is the thickness of the membrane and is constant 


thnoUiSh - out , 


a j b 


U , / If (e^ + + 2 Q ej dx dy (4.6b) 


2 ( 1 ~ V^) 0 0 ^ ^ 


X y- 


4.2.1 Incremental 


j'ls discussed by M-antin (50), in this method, large 
deformations are divided into smaller steps such that during 
each increment of step the problem can be treated as a linear 
problem. Ihe total strain in each step can then be expressed 


(4.7) 

Whereas is the initial strain and|6^| is the developed strain 
during the step. Introducing (4.7) in (4.6b), the strain 


energy becomes 
1 Eh 

^ 2 ^.2. 


/ / [(e°)^ + + zv{e°) (ep + (ep‘ 


( 1 - 0 0 


+ ( 6 ^“=+ 2iJ(e^)Cep + 2(ep(ep + 2(ep(ep 
+ 2v((ep(ep + (e^) (e°^ ta ay { 4 .e) . 




%0 \ 


fhis is further written as 
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v/hcrc 


Uq ^ 


Zs /" •/* ^ ^ ^ (4.9b) 


U 4 =» 
1 


1 Eh 

(1- ■»^) b 0 
a b 


Eh 


(iTv^ { { ®7 ^ ®y 

(4.9c) 


O />0 ^ 8 , / r \ 8 . / i >0 


and 

0, 


Eh 


( 1 - V ) 


- / 2 iy ej e^) dx dy (4.9c) 

2. 0 0 ^ ^ J 


IBae first tern U^ is simply the strain energy present prior 
to imposition of the additional .disturbance' and hence does 
not contribute to the stiffness. 


The second tern IT^ depends on the initial stress and 
hence yields the new stiffness matrix Ik 1 which depends on 

u oj 

the state existing in the elonent prior to the imposition of 
an additional disturbance. This |kgj is known as ‘initial 
stress stiffness’ or ‘geometric* stiffness matrix. 

The third term U 2 depends on the additional . strain. 
This gives the elastic stiffness matrix which is composed of 
the conventional small deflection elastic stiffness T^o] 
the *large displacement* (or initial displacement) matrix 


is a function of nodal disBlacements. In the . 

will not be present . fkl will contain only 


M- M 

present case 

[lcg]ana[kj 

■'■ ■■ 'The' intro;ductioh of r'.'thiS'-';X^ge':,,di^.3.a3C^ent 

in the eq^'libriuii': 
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processes ‘based on trial and error method, for solution. 
Because of dependence of aerodynamic loads on displacements 
and vice versa, this hit and trial method involves lot of 
complexities in numerical computations in a flutter equation, 
where already two unknowns, namely the airspeed and the fre- 
quency, are present in the iteration process. 


iHierefore, as a first step, to simplify the analysis 
yielding a numerically feasible solution, we take 0, 

Purthermore considering the deflections to be ’moderately 
large*, we can derive the, elastic and geometric stiffness , 
m,atrioes in one step onlyt ihus introducing a simple approxi- 
mation. 


Due to the additional in plane extension of the middle 
surface caused by the lateral deflection w of the oscillating 
membrane, the additional strain-displacem’oit relation^ip 

(4.10a) 
(4.10b) 

She additional strains present in ®in, (4.10) represent the 
rotation of the element out of its implane position in the 
xy plane* As stated earlier, the inplane displacements , and 
rotations accompanied With lateral displacements axe not 
considered, ' ■ ■ ■ ■ 


becomes 


a _ /9 Wn' 

'^x *" 2 


and 


*^7 2 '■ij' 
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4.2.2 Derivation ?f pc 


U 


From Eqn. (4.10) and Eqn, (4.9c) v/e pet 
Eh 


^ (1-V^) 0 


a 1; 

/ / 


V-i 

2 

o L 


'^2% ■*■ 2^^^y^x^ 


+ e" (^)‘ 


dx dy 


Noting: that 



E 


i~ y 


1 V 
1) 1 



(4.11) 


Ihe above equation is written as 

u,=lh /Acg(-||)%^;(|E)^a.ay 


(4.12) 


o 0 

writing (4.12) in the following forr., 
a bn„. 1 r o 


u 


1 


1 


/ I 

0 0 


tUL SS. 

?x ay 


0 


T. 


,0 


?w 


KI( dx dy (4.13) 


where and 2?^ are the inplane tensinis in x and y directions 
X y 

respectively, (with ~ ^ hj. s= rr h) 


we have 


'X "" y 
r 


1 a b, , gj _ I 

{ / w [g1 

0 o * ■* 


1 

2 


IE a b $ 

W> C//M 

0 0 


0 


».o 
' T 

y. 

1° o 
o’ 3?° 




ft dx dy 


[G] to ay) {w} 


Gomparing this with ®in, (4.2), we have 

!E 


IT, >= I 



e 


w 


(4.14) 
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where 

^ ^ a h T 

W ■ < i w 


X 


o T' 


0 


[G] 6 .x ay 


(4.15) 


and 



L. 




[G] is 

defined a; 

3 folic 

v:5 (referring to 

Eons, (3.3-3. 5) for 

element 

no , 1 


i X (ct-x) 




L-. 

f,? ® 



W 


ox 


e 


(4.16) 

s: 



i A* (ct-x) 



1 

\ 

Fy 

e 




jii. term hy term integration of Eq.n. (4*15) yields the geo- 
metric matrix pc^j • 

4.2.3 Derivation of |k jj| 


Now third term in the (-.nergy expression (4.9d) becomes 


U. 


Eh 


. / / 

2 ( 1 - 'it ) o 0 


S' -t Am 4 -1 an;. 4 


1 (i|)t (II) dx dy 


‘Px 


ay' 
(4.17} 


Here, because of cubic and fourth order terms gives the 
stiffness matrix as a function of displacements i.e. 

This may be retained or neglected according to the needs of 
the problem. Here it has,.been laei taJcen because of the diffir 
cultios mentioned earlier. So the stiffness matrix pc] for., 
the element will contain [kg] only. Thus for element no. 1, 

M “ [^g] ( 4 * 18 ) 

which is a 12x12 symmetric and complex matrix, fhesss element s^ 
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4 •2.4 Stiffness Matrix for Other SOLeiaents 

Using Eqn, (3.32), the stiffness matrices for other 
elements tore written as follows: 

.5 ^ ..9 . . 13 . . 1 

M * D^J “ W == M given by Eqn. (4.18) 

2 6 10 14 ,• o s.* 1 

W ■= M-M =M = 0-^2 [k]’ 

[k-j = :[.k] = [k] = [k] .= [k] (4.19) 

[k] = [k] = [k] = [k] = e-"6 >-^ [k]’ 


4.3 Derivation of Mass Matrix 

A strai{iit-for;*arcl application of virtual work concept 

gives the mass matrix for element no. 1 as follows* 
i2 ct a b -i2 X x ^ , ,1 


/ / e 
o o 


jfj M ay (4.20) 


M = fn ^ “ 

where is the density of membrane material. It is a symmetric 

complex matrix whoso elements are given in Appendix 1 . 

4.3*1 Mass Matrix for Other Elements ; 

5 9 13 ^ 

£m] « [m] ,= [m] «= [m] given by Eqn, (4*20) 

2 *6 10 14 -i2 'Sa ^^1 

[n] = [m] = [m] = [m] = e ^ [m] 

3 7 11 15 -i4 \a 

[m] = [ml *= [ml - [m] - e , [m] (4.21) 


8 


12 


16 


* 


[m]^ «.{m] M - M -e' 


-i 6 X 3- 


[mj 



OH.-JTER 5 

ASSH4BLY OST EIiSISHTS iiin) ■ffORMULATIOI 
OF EIGBlM¥iJ,Uli PROBLEM 


5*1 Gfcnortil 

Havinr dorive-l -h'’; elemental aerodynamic, stiffness 

and mass matrices, we can put them together to yield the 
equilibrium equation icr an element as 

r 2 *2 2 * 2 ■ 

k - A c n ~ A c A- + i ^ Y q li + V A^ 

: 1 ®2 ..®3. 

^ H (5.1) 

where is a column vector, of generalized forces corres 
ponding t-j the generalized co-ordinates ^wj . Ihe square matrix 
of Eqn. (5*1) represents the generalized ’aerodynamic stiffness 
matrix* for the eloaent. 

Now to obtain a complete solution the two conditions 



(i) displacement compatibility, and 

(ii) equilibrium 


have to be satisfied throu^out. 


(The first condition is satisfied by listing a, system 
of nodal displacements for the whole structure in which 

all the elements participate. Jfext the' conditipn of oyeraU 


equilibrium- within an. el'em.ent is "sal^.ffied 'by ''All 

that is' necessary' ■tO';',sati.efy' 


;60hdi.tion :i4- to', 


li:sh 'eq\iilibrliMi: cGn^i:feione...at'th'e nodes; -of structure* 



4-'6 


*. i. 


Aiiit. iw '".wC . i.j ulic, aoScnblf.;_.;c c.f all the Gloirieirts 

in nliicli the cf^uilibritir! C'ndition of a particulai' no do is 
satisfied by oquaCin/^, the sun of the generalized forces 

contributed by all the clonents neeting at that very node, 
t’. zer", Sius on assonbling the elements by equating the 
correspondin;: renoralized co-ordinates at the global nodal 
pr>intB, a. ’master aerodynaiaio stiffness matrix’ is obtained* 


In essence the assenblat'e of a structure, requires 
the followinr steps: 

1, Transformation of the elemental properties from their 
local Gv'.-ordinatos tc global ones. 


2* Establishment of the transformation matrices W rela- 
tin-f' the elc-nentol no- dal points to the global nodes as 

jvf} = [a] {lf*j (5.2) 

3. Calculation of the characteristics of the structure. 
This is done by assembling according to the congruent 
matrix transformation e,g. for aerodynamic matrices 
wo have 

. [A] = M W (5.3) 

Once the assembly is over the final eigenvalue probloi is 
formulated on app2;ying the homogeneous kinematic boundary , 
conditions by striking out the appropriate rows and columns. 

In the local co-ordinate syst(^ chosen in the present 
case where the' eleoents, lire- rectangular, the el^^jihi'^:;-'', ■■ ;• 
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prr^PGrties in the local raid global co-ordinates system .are 
idcntica,!* The- first stepf therefore} need not be carried 
out here* S ;, we directly move on to the second aspect of 
assdnbly process in the following section. 


5*2 Estn.blishncnt ...f fal Matrices 

In the complete membrane v/e have 16 elements v/ith 
25 global nvdts as numbered in Slgure0.2^ At 'i’th global 
n.)de the generalized displacements are defined as follow's: 


Y/ji^ *= displacement in Z-direction . 


5w 


0 ■= rotation about Y-axis i.e. 

i a-x 


(5.4) 


'5'j, 


rotation about X-axis i.e. 


dw 


and arc- shown in figure 5.1. i takes the value from 1 to 25. 
Thus the (75x1) column vector of generalized displacements 


in global set is defined as 


K! 


w. 


e 


1 


« 

w 




(5*5) 


25 


e 


25 

25 


Similarly for a discrete element the four nodes are npibered 
as shewn in 

3 generaliisod cliap3,ab^'®fe^'^'^' each-nod^. 
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yi I f.?r this cltneiit is defined as 

f vr ] 

! ; I 

f 1 ^’^2 i 


(3,5a) 


■ '^12 


To relate |w| to , we derive the transformation 

matrix *a^j for each element. 


Then for nth element 


[^n] ~ 1,...16) 


i’or elu:jejiit no, I. i 


is i;iivt,n as 


(5.6) 


H 


01 

001 

0001 

00001 

000001 

0. 


001 

0001 

00001 

000001 


ihere 0.,- indicates a row cohtaiJiing 15 coi^s of aeroes 

■. . ; , ' ' . '• .a' '-'r' ’ .i... ■ '« V't' 1' ' .iit* 


and same Is 

other element s . {i# e*‘ ■. ' |^i 6^. ‘ ^i>eh^-x. 
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iCnowin^,? the transfornation natrices for each elenent, 
the .'isser.bly for a particular characteristic (viz. Aerodynamic, 
stiffness, mass) is done according to the equation (5.5) as 
follow si 



16 

£ 

ns1 



1 6 

£ 

npl 


4^n 



(5.8) 


The stiffness matrix [K] thus obtained represents the stiff- 
ness for the entire membrane. 


As is seen from Eqn. (5.7), tho matrices a^*s are 
very sparsely populated containing only unity as the non-zero 
elemented So, while computing the component matrices K^s, to 
store these a^s as such will be a waste of computer space as 
well as of computing time. To avoid this we may ignore all 
zeroes in the specification of the matrix and merely indi- 
cate the position of *undLt’ elements by specifying its row 
number i and column mmber j in the memory location. A typical 
(in Eqn, 5-8) is calculated according to the flow diagram, 
shown in ilg. 5.2, suggested by Argyris (51). 

The first step involves the operation a^^, For a 
specific instruction (1, 3) ^ unit element this requires 

the transfer of the ith elment of k^^ to the Jth column of . 
the space provided for a^» In the second step 

plykjj 

it eleneht ' ■ in question, ' n ^^rahsf er' of the"' 'ith^ row; of "the 
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j th TOW BJii* tllO 0lCD8ll"fc is S*tO2?0d in illG S*toP6 s,llocci*tion 

cl coxed for K# !Phis pi?ocoss is rcpca’tcd for n s= T t;o 16 and 

the clcncnts of ^32*^ arc stored in their respective places# 

In the aaseinbled [k] oatrix a parti oular elaaent will be 

obtained by adding all the elenents in ith row end jth coluian 

of the store allocation contributed by all In this way 

the full 75x75 [Kj matrix is calculated. In a similar manner .• ■ 

the aerodynamic and mass matrices are assembled.^Po achi^e this 
a general assembly scheme is ^own in Hgure 5*44 
5.4 Formulation of Eigenvalue Problem 


While setting up the matrices and during the 
assembly, the membrane is considered to b'e in a free-free 


state in space. For this st at e^ sections 5.2 and 5*3 enable 
us to write the 75x75 master aerodynamic stiffness matrix 


as follows: 


i - A*2 a . ^ 0^ A, + i X* 0 V Aj + TT^ A,] 


(5.9) 


In actual engineering applications such a free-free 
state of meafljrane never, exists- It has to be st^jporfced wi"^ 

some structure in some manner to serve any useful pui^oae* ^e. 
support conditions depend on the needs of the problem- However, 


the master aerodynamic stiffness matrix can be obtained for 
any-' mipport conditions inoa s:^ly -oigitt'ing'-the,, 

appropriate rovs 'is^d ;Odinmns# Is stslJed bi®f@ 2 ?e» -tfee object Pt, 
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flutter, which has bo^.n studied experimentally. 


With thin view, in this analysis the sane confi^ration 
and bruntl'iry conditions are incorporated which were existing 
in the winu tunnel testings* In these experiments the support 
c'?ndition8 had t*.' be chosen as those of fixing the membrane 


at f.ui' if its corners (shown in iigure 5.3) because only 
then it became possible for the flutter to occur within the 
available tunnel speed range. Ihesc conditions correspond to 
aero displaccmonts and rotations at the four nodes namely 
nos* 1, 5, 21 and 25. So, the master aerodynamic stiffness 
matrix for this set of boundary conditions, is obtained by 
striking rov/s end columns number 1, 2, 3, 13, 14, 15, 61, 62, 
63 , 72 , 73 , 74 in the matrix of Eqn. (5.9). 3!he 63x63 matrix 
of Eqn. (5.9) is rewritten as 

K - + M) + i > c V ^2 + -^3 

Iiet (A^ + 


Then above matrix becomes 

ji- ^ i jiI'q y (5,10) 

Generally the aerodynamic dai?>iJig matrix Ig is neglected,, . 
as its contribution is relatively small. But in case It has 
to be retained then we will have to assume the possibility 


that this matrix 


from'll 




multiplier so that we coh write 


w p 


. H-i 
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iVith tl'.i:: v;c ('ut the following equilibrium equation 

|k - ( ^ o'- - ip X c V) + ?^ Aj |f*j «= 0 (5.12) 


Thia c^ui be further written as 


w" [■ 


K + V L 


5 ] [^"1 * ( - ip 


Xc? ) 

(5.13) 


Squo.ti-'n (5.13) represents the conventional eigenvalue 

"^2 2 * 

pr.blen with (Ac - ip Xc ?) as the couplex eigenvalue. 



CHiiPTER 6 

gOLUTIOH AKD QONCLUPri'iirT RT^MARirg 


1 M«:-th., ■■! f Soluti-n 

The clynnnic eiuatiens under the influence of stiffness, 
inertia J'jid -or '.Oynaaic forces have ‘been reduced to the con- 
vtnti-mal f 'rn /f c4f^cnvr.luc problm in section 5.4. Rewriting 
cqunti :n (5*13) ax have 


whfcTe 


and 




jU *s(<iO-ipVCd) 


( 6 , 1 ) 


6^ 2a) 


(6.2b) 


with 6Ja8 the angular frequency which is related to the wave 

# 

velocity c throueii6.'= J^c. (6.2 c) 

The numerical •naly'sio of the systei (6.1) involve following 

steps. 


1. Computation of elemental characteristic ma-iy.ices: 


The aerodynamic matrices can be evaluated after per-^ 
forming the necessary numerical integrations, for the desired, 
panel dimensions *a* and *b*(» The stiffness and mass i^t rices 
are functions of b and A • value of the w^e numb^ X 
depends on the choice' of Wafe lei^h# In the wind tai^el. 
tests, the war®' <^hsa?ved* 


If this 


Information is 
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be . v*ith ?ji wave !!U3do» Kiis then can be 

t.xpj , ■‘.ii a fraction of the stneamwig© disiension of the 
neCibr iiie tiii'.cr tcot# Bins knc-win^? the wave leng’th correspon— 
lir.(i' t ohe iiuoifed panel diiaensions *a* and *b*, the wave 
na-ober A is known. By inserting the corresponding 

vpIU'- f thfoe parameters in the expressions of the aerodyna- 
r.iic, otiffnesa and mass nvatriccs, their numerical values are 
dotcrriinc.'Z • 




bly of I^aJ, [k] r^a [mJ and determination of 


jl-f j and I’K+V^ Aj’j; 

The elemental rar.t rices (k* A and m) obtained in step 
vne are now as&anbled by the method indicated in Section (5,3), 
Then the simple addition of |A^j and [m] gives ^j^Jand that 
of [k] -and gives iij, 

3* Evaluation of matrix [b] : 

Eitrix [b] is obtained by inverting the matrix |l^] 
evaluated in stop two and then p os tmult ip lyin^ +his inverted 
matrix by j^K + 

4# Solution for the eigenvalues of [b]? 

The airspeed T.is the only uziknown contained in 
ooE 5 )lcx matrix |bJ*:,So» once this, is ase^ed> , is ,-k3pw 
completely • men its;;epi^l<». 

by the. standard ' sdbrcft^ihes' o>o;::,o-r,v ; ,,o . 
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Flutt-c r critci i an s 

'She uifcnv'auu, ii, is a function of Y ond^o. The 
('car.il'.x value fi of ^ corr . apendin^^ to a known Y, in general, 
five;: a valu., ?f froQ.ucncy* let the complex eigenvalue 

an,; frtiucricy bo repruoented as 

ju * + i (6.3a) 

and ^ * 6J^ + ihij respectively (6.3b) 

The Teal part ef the frequency represents the oscillations 
while the inaginaiy part (Wj) indicates the amplitude varia- . 
ti'^ns. The condition of instability depends only on i«)j. Whether 
in ncifative or positive, the wave still remains a combi- 
nation cjf sine anf. cosine waves. WhenWj is positive the ampli- 
tude factor becomes thus the amplitude of the motion 

will docruasc with increasing time. When t*)j is negative, the 
opposite is true* 

IfcOj is positive at and negative at Yg (Yg > Y,^), 
then there exists at least one value of Y betwceirv y,j and Yg, 
at whidaiOj vanishes* At this speed is real, corresponding 
physically to a simple harmonic motion, Srch a speed will 
separate the speed range in its nei^borhood into two regions, 
in one of which^I > ^ where the motion is damped and stable; 
in the' otheriO^j ^ 0 whe?e the -aaplitirde, increases, . 
and the motion ,is unstable* 
called t,he . * critical ' *■ ' :■ ' : r : i ■■■ 
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S r-ii- iv. ccncerncd, tw;. cases are possible 

it th. .;ritica,l n Either vanishes, or it does not 
V ** “hvn the displacements are independent of 

tl::;,., : ut tra. otructuro has lost its power to recover its 
■ritin-d 1 .i..j **hcn oisturbed* The panel is said to be in cri- 
tic-a divcr^-Ult coniLtion. 0, the motion is harmonic 

v.'ith a.n iiiiiiil init<. nraplitudc. It is ssiid to be in the critical 
rlutt.^r c -nlition, Hcncc, 

(Oj «= 0, IS 0 Implies divergence 

« 0, ^ 0 implies flutter. 

In both cases the aerti clastic system nay be said to be neutr- 


ally stable. 


Prom the available experiments, it looks as thou^ 
has a very high positive value initially# As the speed is 
raised gradually, it keeps on decreasing till the onset speed . 
is reached, where it becomes zero. At a aligh'tly hi^er speed 
oscillations become divergent and it seemsiSJ^ becomes negative 
iiD the phenomenon of divergence was not observed, it is felt 
therefore thatiilg^ a non-zero value at the critical speed 
and the instability is recognised as the flutter. 


as 


Now using lan# (6#3) the B(ln. C6.2b) is rewritten 

' '■■■ ' ' \ ■ :Z- ■■ '■ " " 


Applying tah® :abor‘t!f%tte3f/'o^ltarionr,^ 
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% ^ ^ P ^cr 


F/luitinr thi. rvnl nnd inp.rinnry partn,- we get 

1/2 

1/2 

.Jirl >»i- “P^cr ^R> 


(6.5) 


( 6 . 6 ) 

(6.7) 


Bqn. (6.7 ) {■ivcss the relsition between liie reo^l and iniaginaj’y 
parts :>f the eigenvalue » at the critical condition. 


The values of airspeed Y is now gradually increased 
and the corresponding (lowest) eigenvalue is found. The process 
is Continued till the eigenvalue obtained is such that it sati- 
sfies the flutter criterion (6.7). The airspeed (7 ) and 

the frequency (^3.®= corresponding to this condition 

give the flutter speed and the flutter frequency respectively. 
If the aerodynamic damping is. neglected, then Eqn, (6.2b) redu- 
ces to 




2 

60 


and the corresponding flutter criterion becomes 




A Lcs 

i^rZ- 


and 


IIt * 


•fw 


( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 


Hence here we will seelCA'tii® existence of real eigenvalue for 


the flutter to occur* 

The method of solution described above in stops 1*^4 •: 

can be. repi'esenfced by a tOTical flow chart as sho^' below. , 























anrilyticrCL foiniulation for noderately large anpli- 
tu<lc tr.avcllijig wave nonlrane flutter in inconpressible flow 
h":3 been prcBcnted using finite olenent technique. The flow 
hao been assuned to take place on both upper and lower sur- 
faces of the racDbrone having flexible bound:iries. 

Unfortunately the aforementioned problem is of consi- 
derable complexity because of high order . complex matrices 
encountered in the finite element method, structural nonlin- 
u-.'ixitic'S and complicated aerodynamic expressions. The numerical 
evaluation of the integrations involved in the establi^ment 
' f aurodynanic matrices and the computation of eigenvalues of 
the cemplc.x matrix is undoubtedly a difficult and time consu- 
minr task* In the present thesis, only the analytical forrau- 
lati »n could bo dealt with. Therefore, the investigation is 
by n.) means a ctxipleto one without the numerical analysis. 
However, in the formlation, every possible attempt has been 
made to present it in a form which is readily adaptable to , 
computer programming. This helps in reducing the bime .required 
durinf* the ijrulininrry adjustments of the computer prograraaing. 
The formulation presented is made general as for as possible, 
to account for any set of boundary conditions, for a rectan- 
gular membrane. The specific boundary conditions chosen cor- 
respond to fixing the mmibTme at its four corners, as thts 
is the situation existing in the experimental study of the 
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I r .i.lcij* Thlo 't 5 rpiGfil coiii'ipuii'a'ti,;)!! 00111101; Ido analyssd liy 
r.n^ a38Une.I m \c solution, ac no sinple function is able to 
satiijfy these boundary conditions* !Ehis reveals the necessity 
)f finite element method because of its capability to treat 
such priblcns* 

The naabrano ho.s been divided into 16 rectan^lar 
eloncnts fornihj; a 4x4 Cfidwork, A 12 parameter non-conforming 
model has been used* Small perturbation theory for non-circu- 
latory inoonprcssible flov/ has been employed to give the aero- 
dynamic forces. The nonlinearity caused by the large deflections 
has been t'^Kon cai*e of by deriving the geometric stiffness 
matrix. The aerodynamic matrix involves the integrations which 
h-'.v«; tw be evaluated numerically. The elements for stiffness 
and man a matrices have been found by exact integrations. Their 
numerical values can be obtained on direct substitution of the 
panel dimensiuna and vvavo number. 

The effect of inplane tension, size, aspect ratio and 
material density, on flutter behaviour can be readily studied 
by varying the appropriate parameters without altering the entire 
formulation. 

6.3 Further Scope 

To complete the investigation the first obvious ta^ 
will be to find the numerical solutions for the present foriau- 
latlott. Once the results are obtatoed they can be ccmip^ed wl# 
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thw uxi-criDcntn,! rusults and the necessary inprovenents in 
tirsr analysis t:rin be Icne, if needed, 

Hrv/cver, the follovi;ing extensions are worth carrying 
.;ut tu strengthen the analysisi- 

(1) Incorj.orating the cenforning nodel in the analysis* 

(2) Increasing the nunher of elements. 

( 3 ) Applying the incremental step approach with more 
launher of steps by reducing the step size. 

(4) Considering the large displacement matrix • 

(5) Censidering: the cxirved shell elements during process 

( 2 ), 

(6) Modifying the aerodynamic theory to study the, flutter 
bchavi .'ur of the membrane placed at an angle of attack » 
and at hifher speeds. 

iUiy other configuration can be analysed by using the 
triangular elements, piong the some lines as discussed in 
the present work. 
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